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The electric dipole polarizabilities and hyperpolarizabilities for the lithium isotopes 
the ground state 2 2 S and the excited states 2 2 P and 3 2 D, as well as the leading resonance and dis- 
persion long-range coefficients for the Li(2 2 S)-Li(2 2 S) and Li(2 2 S)-Li(2 2 P) systems, are calculated 
nonrelativistically using variational wave functions in Hylleraas basis sets. Comparisons are made 
with published results, where available. We find that the value of the second hyperpolarizability of 
the 2 2 S state is sensitive to the isotopic mass due to a near cancellation between two terms. For 
the 3 2 D state polarizability tensor the calculated components disagree with those measured in the 
' sole experiment and with those calculated semi-empirically. 
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I. INTRODUCTION 



The energies, transition probabilities, and polarizabilities of the lithium atom and its isotopes and their mutual long- 
range interaction coefficients have been extensively studied, and many calculational approaches have been developed 
£>y and tested against each other and against available experimental results. Nevertheless — as we shall show — there is 
i-£h ' still a need for calculations of increasingly high-precision to serve as benchmarks, to predict atomic and interatomic 
coefficients, and to help understand discrepancies between various results. Calculations at the level reflecting the 
mass of lithium isotopes may be applied to a diverse set of recent areas of interest in, for example, astrophysics 1], 
t— I • ultracold atom-atom scattering @, y] and Feshbach resonance analyses [1, photoassociation spectroscopy [f| 7], 
J> ' atom-molecule scattering Q and three-atom inelastic collisional loss studies Measurements and calculations on 
excited states and their properties are also at the frontier, though discrepancies between theory and experiment remain 
even for the lowest fine structure levels [ToL ITlT ] . 

It has been demonstrated that ab initio wave functions obtained variationally using Hylleraas-type basis functions 
are capable of yielding highly accurate results for Li properties, cf. [l2l . [l3j . In turn, those results have proven useful in 
gauging the effectiveness of Gaussian-type basis functions [3| and other calculations [H| , semi-empirical methods [l(| , 
and relativistic approaches [U Ell flil. The purpose of the present work is to apply and extend methods developed 
over a series of previous papers (e.g. [I2l.l20l.l2ll]) to the excited 2 2 P and 3 2 D states of the lithium isotopes thereby 
providing a consistent and highly accurate ab initio treatment of the polarizabilities and their related quantities using 
the nonrelativistic Schrodinger equation. We also refine the previous results [l2j for the 2 2 S ground state by improving 
the accuracy and by the inclusion of the isotope mass. For the excited states of the isotopes we provide calculations 
of static polarizabilities and static second hyperpolarizabilities and we investigate the excited state resonance and 
dispersion long-range interaction potential energy coefficients. 



II. THEORY 



In this section, the Hamiltonian and basis sets will be exhibited along with expressions for the electric multipolc 
transition operators, polarizabilities, second hyperpolarizabilities, and dispersion coefficients. The formulation extends 
and unifies those given previously for Li [l2l.l20l.l2ll]. for He [22|, and for (23|; we include the isotopic mass and 
we treat the response of the atom to an applied electric field. 
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A. Hamiltonian 



The transformation from the laboratory frame to the center of mass frame that we will use for describing the Li 
atom was given by Zhang and Yan [23| in a general form for n + 1 charged particles. It was applied to generate the 
Hamiltonian and transition operators for H^ in Ref. [23| and for He in Ref. [22j . 

We directly follow the expressions given in Eqs. (20)-(26) of Ref. (23|, where the particles are labeled by the index 
i = 0, ...,71. We make the replacement n — 3, identify particle with the nucleus of mass m , identify the particles 
i = 1, 2 and 3 with the electrons, so that mi — m 2 = 777,3 = m i where m is the electron mass, to obtain the Hamiltonian 
expressed in the center of mass frame describing the internal motion of the Li atom, 



i=i *>j>i i =i *>j>i 

where is the position vector of electron i from the nucleus, and fi = mmo/(m + mo) is the reduced mass. In this 
subsection j stands for a summation index (in subsequent parts of the paper it will be an angular momentum quantum 
number), and qj, j — 0, 3 are the charges of the respective particles. 
The 2 £ -pole transition operator is 

3 

T e = Y,H R i Y m{%) , (2) 

i=0 

where R;, i = 0, ...,3, as shown in Eg. (27) of Ref. [23j. It can be transformed into the center of mass frame by 
applying Eqs. (23) and (24) of Ref. [23|], which may be written in the form 

3 

with €ij = Sij — mj/Mr, i = 0, 1, 2, 3, j — 1, 2, 3, and Mr = mo + 3m. The general formula for Ti expressed in the 
center of mass frame is given in Eq. (30) of Ref. [23| . 

For a four-particle system, it can be shown that the transition operators Tg with t up to 3 have the following explicit 
forms: 

Ti = feft^Wiofo). (4) 
j=l ^ i=0 ' 

T2 = E (j^^r^oiv,) + (Eftw)^*^.®^ , (5) 

j=l \ i=0 ' {j,fe} \ i=0 ' 

with {j, k} = (1,2), (2, 3), and (3, 1), and 

T 3 - E(E«4)r? y »(*i)+sJ§ £ (E»^)^((*i®'i) (a> ®^ 8) 

j=l ^ i=0 ' j,fc=l *=o ' 

^ 4=0 ' 

For a neutral system, the finite mass effect enters into the transition operator Tg as a polynomial of degree t — 1 in 
rrij/MT- For Li, the three coefficients appearing in the sets of parentheses in Eq. ^ for the transition operator T3, 
for example, can be written as 

E«4-«-^|£) + M3:V- CO 



i=0 

3 



;=o 
3 



i=0 
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For an atomic system of infinite nuclear mass, all the coefficients are zero except for those of rjY^Tj), which are 
equal to qj, as expected. We now let m = 1 to utilize atomic units for the remainder of the paper. 

B. Variational basis sets 

The wave functions are obtained by solving variationally the energy eigenvalue equation for lithium 

ff *o(ri,r 2 ,r 3 ) = £ *o(ri, r 2 , r 3 ) , (10) 
with Ho given by Eq. ([I]), in terms of the explicitly correlated basis functions in Hylleraas coordinates, 

^ra.ra) = rf rf rf rj^r^^e-^-^-^yJ^^^.fa^sMl.Z.S) , (11) 

where = |r, — r ^ | is the inter-electronic separation, y^l 2 y 12 1 3 is a vector-coupled product of spherical harmonics 
to form an eigenstate of total angular momentum L and component Ml , which can be written in the form 



all mi 

x (£ 12 mi 2 ; £3m 3 \e 12 e 3 ; LM L )Y eimi (ri)Y e2m2 (r 2 )Yi, 3m3 (r 3 ) , (12) 

and x(l,2,3) is the three-electron spin 1/2 function. The variational wave function is a linear combination of 
basis functions 4>, anti-symmetrized. The procedures followed are similar to those described in Ref. [2(|. With some 
truncations to avoid potential numerical linear dependence, all terms in Eq. (jlip are included such that 

3i + h + h + jl2 + j 2 3 + hi < ^ . (13) 

where £1 is an integer, and the convergence for the energy eigenvalue is studied by increasing progressively. The 
basic type of integral that appears in this work is of the form 

dndr 2 dr 3 rf rf rf rf 2 2 r 3 2 f 4i e ~ ari _/3r2 " 7r3 
x Y e \ m , (h)Y e l m , (r 2 )Y e \ m , (h)Y hmi (h)Y iama {h)Y tams (h) (14) 
and computational details for this integral can be found in Ref. [2lj |. 

C. Stark effect and polarizabilities 

The polarizability of an atom can be considered as a measure of the response of the charge cloud to an external 
electric field, which can be illustrated by the Stark effect. Consider a lithium atom in a weak external electric field 
£ = £ i. The initial state is assumed to be a parity eigenstate and is written as |0) = \uqLM), where no is the 
principal quantum number and L and M are the usual angular momentum quantum numbers. According to the 
perturbation theory, the energy shift due to £ can be expressed in the form 

AE = AE 2 + AE 4 , (15) 

where AE 2 and AE4 are, respectively, from the second- and fourth-order corrections, whereas the first- and third- 
order corrections are zero because of the parity selection rule. The detailed derivation for AE 2 and AE4 is given in 
the Appendix. The final expression for AE 2 is 



£ 2 

AE 2 = 4 



(16) 



ax+a { pg 2 (L,M) , 

where g 2 (L, M) is the only M-dependent part, defined in the Appendix by Eq. (|A29|) . and ot\ and are, respectively, 

(T) 

the scalar and tensor dipole polarizabilities. The polarizabilities at\ and a\ can be expressed further in terms of the 
reduced matrix elements of the dipole transition operator: 

cm = y^ai(Xq) , (17) 
a { P = Y,W(L,L a ) ai (L a ), (18) 
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where 



= 8^ v ^LWT.WnL^ 
1[ al 9(2L + 1) ^ E n (L a ) - E no (L) 



(19) 



with T a = ^^ 0% i? 4 y 10 (R, ; ), Eq. ©, and 



(2L + 3)(L + 1) 



i L L, 



(20) 



In the above, the set of energies and wave functions {E n (L a ), \nL a M a }} correspond to an intermediate energy spectrum 
allowed by the dipole selection rule, which can be obtained by diagonalizing the Hamiltonian in a Hylleraas basis set 
of given symmetry L a . In particular, for the case of L = 0, 



= oi(P) 

, = 0; 



Oil 



for L = 1, 



ai(S) + ai(P) + a^D) , 
-a 1 (S) + ~a i (P)-±a 1 (D) 



and for L 



Oil 

AT) _ 



ax(P) + ax(D) + ai (F) , 
-on(P) + ai(D) - ~a x {F). 



(21) 
(22) 

(23) 
(24) 

(25) 
(26) 



In Eqs. (|23|) and ([24|l . ai(-P) is the contribution from the even-parity configuration (pp')P. In Eqs. ((25|) and (|26| . 
ai(-D) is from the odd-parity configuration (pd)D. 

The fourth-order energy shift can be written in the form 



£ 4 

AE 4 = - — 
24 



7o + 72 52 (£, M) + 74 34 (i, M) 



(27) 



where g<t(L : M) is given by Eq. (|A47[) in the Appendix. In Eq. (|27[) . 70 is the scalar second hyperpolarizability, and 
72 and 74 are the tensor second hyperpolarizabilities, which can be written as 



To 



(-1) 



2L 



128tH 



/ ^ Go{L,L a ,L b ,L c )T(L a ,L b ,L c ), 

V 2L + 1 



(28) 



LaLfrLc 



72 = (-1 



,2i 



128tH 



L(2L - 1) 



(2L + 3)(L+ 1)(2L- 



La.Lf,L c 



(29) 



74 = (-1) 



2L 



3 V (2L + 5)(L + 2)(2L + 3)(L + 1)(2L + 1) , 4< 



L a L b L a 



where T(L a , L bl L c ) and Ga{L 1 L a , L b: L c ), respectively, are defined in the Appendix by Eqs. (|A44|) and (|A51[) . In 
particular, for the case L = we only need to consider 



7o 



1287T 2 



-|-T(1,0,1) + ^T(1,2,1) 



(31) 
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[Note that the case Lf, = 1 does not enter in Eq. (|3Tj) . The first 3-j symbol in (|A51[) requires K\ — or 2 and since 
L = here, the first 6-j symbol requires = K\] therefore, — or 2.] For the case L = 1, we have 



7o 



1287T 2 



3 
^5 



^T(0, 1, 0) + Ar ( 0, 1, 2) + lr(l, 1, 1) - lT(l, 1, 2) - i_T(l, 2, 1) 



72 



45d T(1 < 2 ' 2) + il r(2 < x ' 0) - ^ T(2 ' x) + iio r(2 ' 2) - i T(2 < 2 ' x) 

^r(2,2,2) + A r(2 ,3,2) 

128712 [ - ^ °) - il^ ' 2 > + is 7 * 1 ' m) - ^a, i, *) - so^' 2 ' 1} 

— T(l, 2, 2) - — T(2, 1, 0) - — T(2, 1, 1) — T(2, 1, 2) - — T(2, 2, 1) 

900 v ' ' ' 135 v ' ' ' 180 v ' ' ' 2700 v ' ' ' 900 v ' ' ' 

— T(2, 2, 2) — T(2. 3, 2) 

900 v ' ' ' 1575 v ' ' ; 



(32) 



(33) 



and for L — 2, we have 
128tt 2 1 



7o 



140 T(l, 0, 1) + 119 T(l, 2, 1) + 84T(1, 2, 3) + 84T(3, 2, 1) + 74T(3, 2, 3) + 60T(3, 4, 3) 



3 15750 

105 T(l, 1,1) + 2lV5T(l, 1, 2) + 21 T(l, 2, 2) + 2lV5 T(2, 1, 1) + 21 T(2, 1, 2) + 21 T(2, 2, 1) - 119 T(2, 2, 2) 



56 T (2, 2, 3) + 56 T(2, 3, 2) + 4V70 T(2, 3, 3) + 56 T(3, 2,2) + 4V70 T(3, 3, 2) + 20 T(3, 3, 3) 



(34) 



1 



72 



1960 T(l, 0, 1) + 1225 T(l, 2, 1) + 840 T(l, 2, 3) + 840 T(3, 2, 1) + 380T(3, 2, 3) + 240 T(3, 4, 3) 



1287T 2 

3 ' 154350 

735 T(l, 1, 1) + 147\/5T(1, 1, 2) + 147 T(l, 2, 2) + 147V5T(2, 1, 1) + 147T(2, 1, 2) + 147 T(2, 2, 1) + 1519T(2, 2, 2) 
448 T(2, 2, 3) - 448 T(2, 3, 2) - 32 V / 70T(2, 3, 3) - 448 T(3, 2, 2) - 32\/70T(3, 3, 2) - 160 T(3, 3,3) 1 , (35) 



1287T 2 



1 



74 = 



490T(1,0,1) + 49T(1,2,1) + 14T(1,2,3) + 14T(3, 2, 1) + 4T(3, 2, 3) + -T(3,4,3) 



3 128625 

245 T(l, 1, 1) + 49V5T(1, 1, 2) + 49T(1, 2, 2) + 49%/5 T(2, 1, 1) + 49 T(2, 1, 2) + 49 T(2, 2, 1) + 49T(2, 2, 2) 



14T(2, 2, 3) + 14T(2, 3, 2) + V70T(2, 3, 3) + 14T(3, 2, 2) + V70T(3, 3, 2) + 5 T(3, 3, 3) 



(36) 



In each of Eqs. (|34|) - (|36[) . the terms in the first set of square brackets only involve the intermediate states of natural 
parities, which make the dominant contributions to the hyperpolarizabilities, while the terms in the second set of square 
brackets involve the intermediate states of unnatural parities, which make subordinate contributions. For example, 
the term T(1,0, 1) involves the intermediate states ls 2 np 2 P°, ls 2 ns 2 S e , and ls 2 np 2 P°, which are all natural-parity 
states, while the term T(l,l,l) involves the electronic configurations ls 2 np 2 P°, lsnpn'p 2 P e , and ls 2 np 2 P°, where 
two states are natural-parity states and one unnatural-parity. 

The scalar dipole polarizability defined in Eq. (fT7|) can be generalized to the 2 -pole polarizability at 



(37) 



where 



(T s = 8^ v \(n L\\T e \\nL a )\ 2 

ae[ a) (2£+l)a(2L+l)4- E n (L a ) - E no (L) ' 



(38) 



and Ti is the 2 f -pole transition operator given by Eq. @. 
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D. Coefficients for long-range interactions between two atoms 

First, let us consider the simplest case where both a and b are Li atoms in their ground states [12]. At large 
separations R, using second-order perturbation theory, V ab can be expressed as a series in inverse powers of R, 

V ,-_9l-9l-9^. (39) 

Vab ~ i?6 RS RIO ' ^ 



where 



C 6 = -G ab {l, 1) , (40) 

7T 



and 



C 8 = ^-G ab (l,2) + ^-G ab (2,l), (41) 



C 10 = —G ab (l, 3) + — G Qb (3, 1) + — G afc (2, 2) . (42) 

7T 7T 7T 



Introducing the oscillator strength for the transition \uqL) — > |nL') 



^ = (2£+l)^2£ + l) £; " o|(noL||rf||n£,)|2 ' (43) 



where £"„o = E n (L') — E no (L) is the corresponding transition energy, G ab (i a ,i b ) can be written in the form 

fit a) 7(h) 

nO Jn'O 
2 frf E nO E n 'o( E nO + E n'o) ' 



7T / / / 



where throughout a and 6, respectively, denote atom a and atom b. 

Next, we consider two like lithium atoms a and b, where atom a is in the ground state and atom b in an excited state 
with orbital angular momentum L b and associated magnetic quantum number M b . The zeroth order wave function 
for the combined system ab can be written in the form [l2l |: 

* ( °) - ^[^ a (a)^ b {L b M b -p)+l3^ a {p)^ b {L b M b -a)] , (45) 
V2 

where er and p represent, respectively, the set of all of the internal coordinates for atom a and atom b, and /3 = ±1 
describes the symmetry of the system due to the exchange of two atoms. According to the perturbation theory, the 
first-order interaction energy is given by 

y( 1 )(L b M fe ;^)=-^±i, (46) 

where 

<&l = P(~l) 1+LM (2L ^ 1)2 ( Lb 2 + b Mb ) |(*aH||T ib ( CT )||* 6 (L 6 ; CT ))| 2 . (47) 

One can see from Eq. ([46]) that, for the Li(5)-Li(P) system, the interaction energy is proportional to R~ 3 . To get the 
next order energy, let the complete set for the intermediate states of the system be 

{Xs{L s M s -cr)u t {L t M t -p)}, (48) 

with the energy eigenvalue E^' = E^ + E^ . According to the second-order perturbation theory, the second-order 
energy is 

V R? RS { yj 
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where 

C t = E . (5°) 

= y n*. — (51) 

and the energy for the unperturbed system is E^ = E^ + eP\ Following one can obtain the following 
expressions for fig* and fig* that are in agreement with the formulas in Ref. [221 ]: 

nf = |<* ( t r)||T 1 («r)||x.(l; a)) | 2 £ G(l, 1, 1, A, 1, M 6 )|(* 6 (l; p)||Ti(p)||o; t (A; p))\ 2 , (52) 

A 

ilf = 2 |(^ tt (^)||r 1 (^)|| Xs (l; C r))| 2 ^C(l,3, 1, A, 1, 7V4r b )<^ b (l; J p)||r 1 ( >0 )||^ f (A; / >)> 

A 

x<* 6 (l;p)||T a (p)||w t (A;p)) 
+ |(* ( -)||T 1 ( -)|| Xs (l;<T))| 2 ^G(2 ) 2a,Aa,M 6 )K* 6 (l;p)||T 2 (p)||^(A;p))| 2 

A 

+ |(* ( -)||T 2 ( -)|| Xs (2; CT ))| 2 ^G(l,l,2,Aa,M 6 )K* 6 (l;p)||T 1 (p)||^(A;p))| 2 

A 

+ G3(2,2,l,l,l,M 6 )(* a (<r)||ri(<r)||x.(l;<r))<* 8 (p)||Ti(p)||w t (l;p)) 

x<* 6 (l;p)||T a (p)||w t (l;p))<* 6 (l;o-)||T a ( < r)||x.(l;o-)) 
+ 2G3(l,2,l,2,l,M 6 )(* (<r)||ri(o-)||x.(l;<r))<* a (p)||r 2 (p)||w t (2;p)) 

x (* 6 (1; p)||T 1 (p)||w t (2; p))(* 6 (l; cr)||T 2 (fT))|| Xs (l; <r)> 
+ G 3 (l, 1,2,2, l,M 6 )(* (<r)||T 2 (<r))||Xa(2;ir))<* (p)||T 2 (p)||« t (2;p)) 

x (* 6 (1; p)||T 1 (p)||w t (2; p))(* 6 (l; -)||T 1 ( -)|| Xs (2; <r)) , (53) 



with 



G(L,L',L s ,L t ,L b ,M b ) = ^ (L,L'y 1/2 E 

1 S + > Ms M t 

Lb L Lt \ l Lb L' Lt 
-M b M s M t I \ ~M b M s M t 



and 



G,{L,L',L S) L u Lb,M b ) = (-l) L+L ' j- (L,L')- 1/2 E (~1) M ' +Mt ^ K^ t 

(2L. + l)(2£ t + l) ^ 

Lf, L L t \ I Lb V L s 
-M b M s M t ) \ -M b M t M s 

In Eqs. (|5ll) and ([SB"]), the coefficient is 

-,1/2 











A'- 






t" p A 



and £,•■■) = (2£+l)(2L + l) 



III. RESULTS AND DISCUSSIONS 



(54) 



(55) 



(56) 



In this section, we present the results of the calculations of the static electric dipole, quadrupole, and octupolc 
polarizabilities, the second hyperpolarizabilities, and the resonance and dispersion long-range coefficients for the pairs 
of atoms using the wave functions obtained as described in Sec. Ill Bl In addition, we give some calculated oscillator 
strengths and sum rules that might be useful. 
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A. 2 2 S state: Polarizabilities and hyperpolarizabilities of atoms and dispersion coefficients between two 

atoms 

TableUpresents a convergence study for the calculations of the scalar dipole polarizability ot\ of lithium with infinite 
nuclear mass, °°Li, in the ground state. In the table, Nq and Np, respectively, are the sizes of the basis sets for the 
ground state and for the intermediate states of symmetry P. The extrapolation, obtained by assuming that the ratio 
between two successive differences in a\ stays constant as the sizes of the basis sets become infinitely large, yields 
the value 164.112(1). This is in perfect agreement with the value 164.111(2) of Ref. [l2j, based on calculations up 
to much smaller values Nq = 919 and Np — 1846, confirming the efficacy of the extrapolation method used in that 
work. For 6 Li and 7 Li, a similar convergence pattern exists. 

Calculations on the hyperpolarizability 70, on the other hand, require much larger basis set sizes to achieve accuracies 
even approaching that achieved for the polarizabilities. In our approach, there is a partial cancellation of significant 
figures between the sum of the two terms | T(l, 0, 1) and 2, 1) in Eq. (f3"Tj) . even though the individual terms 

are converged to about 4 significant figures. Table HT1 presents the convergence study for calculations of the hyperpo- 
larizability 70 of °°Li in the ground state. At the the largest sizes of basis sets in Table 02 §T(1,0, 1) = -3463.861 
and T(l, 2, 1) = 3471.078, resulting in a loss of about two significant figures when added. 

Table [nil summarizes the final values of the polarizabilities ot\, a.2, and 03, 70, and the dispersion coefficients C§, 
Cg, and C10 for °°Li, 7 Li, and 6 Li in their ground states 2 2 S. The dispersion coefficients were calculated similarly to 
those calculated in Ref. [12l ]. but with the treatment of finite nuclear mass. In order to maintain numerical stability, 
especially in diagonalizing the Hamiltonian for large basis sets, all the calculations were performed using the multiple 
precision arithmetic QD developed by Bailey et al. (24j, which has 64 decimal digits. The finite nuclear mass effect is 
most significant for 70, as shown in Table 11111 Fitting the results, we find, roughly, that 70 ~ 3060[1 — 1000(l/mo)], 
where mo is the nuclear mass and where the coefficient of 1/mo is about 1000 times larger than it is for, say, ot\. The 
effect is numerically significant since, as discussed above, our convergence studies show that at least two digits of the 
hyperpolarizabilities are accurate. It would be interesting to explore this effect for other systems. 

In Table HVl we compare our results to some of those from the literature for the polarizabilities a±, a.2, and a?, 
and for the dispersion coefficients C@, Cs, and C10 of ground state Li atoms. As most of the previously published 
works have been devoted to °°Li, we list our infinite nuclear mass calculations in Table IIVI The most accurate 
Hylleraas-type calculation prior to the present work is that of Ref. (l2j |. which was obtained using a basis size up 
to 1846 using methods similar to the present work, extrapolated as discussed previously. An extensive tabulation of 
over forty results for a\ covering much of the published work from 1959 to 1996 can be found in the review article 
by King [l3T ] and another tabulation is given in Ref. [12l| . Tabulations including ot2 and 03 are given in Refs. (l2| 
and [25[. In Table HVl we tabulate some theoretical results from between 1996 and the present. The exponentially 
correlated Gaussian-basis set calculations of Komasa [26j for a\ are in excellent agreement with the present work and 
they were carried out with a much smaller basis size. There is a slight discrepancy, however, for a2, but Ref. [2fj| does 
not provide computational uncertainties. Both a\ and «2 were calculated using semi-empirical model potential-based 
methods by Cohen and Themelis (2?| and by Zhang et al. [liij]. Cohen and Themelis [27| used a method dubbed 
Rydberg-Klein-Rees quantum defect theory (RKR-QDT), which utilized experimental energy levels and it contained 
some adjustable parameters fixed using the experimentally determined 2p state lifetime. The results of Zhang et 
al. [Tfij were obtained in the framework of a frozen core Hamiltonian with a semi-empirical polarization potential. 
Of the two semi-empirical approaches, the RKR-QDT results are in much better agreement with the present work; 
the results of Zhang et al. for 0:3 vary from the present work by the same percentage as do their calculations of ot\ 
and ci2- Chen and Wang [25| evaluated 02 and a-$ for the ground states of lithium-like ions using the full core plus 
correlation method. The present values lie just outside the lower limit of Chen and Wang's error bars, but Ref. (2f| 
does not reveal how the uncertainties were obtained. 

We can also compare our results to recent relativistic calculations. For a\ generally the effect is to reduce the non- 
relativistic value by a term of order (1/137. 037) 2 . Derevianko et al. 28] and Porsev and Derevianko [29| calculated a±, 
a.2, and 03, as well as Cg, Cs, and C10 for the ground state of lithium using relativistic many-body perturbation theory 
(MBPT). Sahoo et al. [lalSO] performed ab initio relativistic coupled-cluster calculations on the dipole and quadrupole 
polarizabilities, and the dispersion coefficients Cq and Cs of Li. Very recently, Johnson et al. [rp also calculated 
the polarizabilities of 7 Li by applying relativistic MPBT. For ax, the results of Refs. [28[ and Il7ll . obtained using 
relativistic MPBT are smaller than our results by the expected factor. The results of Refs. [HI, [3(| are significantly 
lower, perhaps due to neglected higher-order effects [l8[, see Table ITVl 

For Cq, compared to the previous value the uncertainty in the present result has been reduced by a factor of 3. 
The semi-empirical calculations of Zhang et al. and the relativistic calculations of Refs. [28| and [2^] are in good 
agreement with the present results, though Zhang et al. slightly overestimate the coefficients, as shown in Table HVl 

The most precise measurement of a\ was that of Miffre et al. f3l"j obtained with 0.66% uncertainty using atomic 
interferometry. It is a factor of three more precise than the earlier measurement of Molof et al. [32| • These values are 
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included in Table EH 

Currently, however, the experimental accuracy has not reached the stage where finite nuclear mass and relativistic 
effects in ai can be tested stringently. One can see from Table HTT1 that the finite nuclear mass correction to a\ for 7 Li 
is 0.049 a.u., of which 0.038 a.u. comes simply from the mass scaling of the Bohr radius. The remaining part 0.011 a.u. 
is due to the mass polarization terms (— 1/mo) Vj • Vj in the Hamiltonian Eq. |T]), where mo is the mass of the 7 Li 
nucleus. Furthermore, although the relativistic effect has not been evaluated rigorously, it can be, however, estimated 
to be —0.06 a.u. based on a relativistic coupled-cluster approach [U [34], [35J. As for the QED effect, Pachucki and 
Sapirstein [36| performed a relativistic and QED calculation on the dipole polarizability of helium and found that 
the QED correction is about a factor of 2.5 smaller than the relativistic correction and is opposite in sign. If we take 
this reduction factor for the case of lithium, the QED correction is thus estimated to be 0.02 a.u. Hence, the finite 
nuclear mass, relativistic, and QED corrections cancel out almost entirely, just as in the case of helium. However, 
a definitive conclusion regarding this cancellation can not be drawn until the calculation of Pachucki and Sapirstein 
can be extended to the case of Li. 

As discussed above, we have found that the hyperpolarizability is extraordinarily sensitive to the finite nuclear 
mass. In Table [V] we compare our calculated value of 70 with some of the published results, all for 00 Li. A more 
comprehensive table of earlier work is given in Ref. (l3j . Pipin and Bishop [37| calculated 70 by applying the combined 
configuration-interaction-Hyllcraas method. Their result, with one significant figure, is in good agreement with ours. 
Kassimi and Thakkar 38, 39] used the coupled-cluster approach, where the reported uncertainty of 10% in 70 was 
later re-adjusted to 3% [39l |. Laughlin (40| performed a semi-empirical one-electron model potential calculation and 
found that the final result for 70 was highly sensitive to the data used, particularly the value of ct\. The value 
from Jaszuhski and Rizzo [4l[ was obtained using a series of multi-configuration SCF wave functions. Finally, Cohen 
and Themelis [27| computed 70 using RKR-QDT. Comparing to our result, the level of accuracy they achieved is 
about 10%. The RKR-QDT calculation was sensitive to the potential adopted, as discussed in Ref. [27|- The present 
methodology has the advantage that no adjustment is required. Once the convergence pattern is established, the 
extrapolated value should be reliable. 



B. 2 2 P state: Polarizabilities and hyperpolarizabilities 

Table IVTl shows the convergence of ct\ and of calculated for °°Li in the 2 2 P state, where N$, Nrp p t\ P , and Nd 
are the sizes of bases for the intermediate states of symmetries S, P, and D respectively, and (pp')P stands for the 
main configuration of two p-electrons coupled to form a total angular momentum of P. Since the contribution from 
the (pp')P configuration is well converged at N^ pp ^p = 3413, we did not increase N(ppi\p any further. 

Table IVHl summarizes the final values of the scalar polarizabilities a%, cx.21 and 03, the tensor dipole polarizability 
the scalar hyperpolarizability 70, and the tensor hyperpolarizability 72 for °°Li, 7 Li, and 6 Li in their 2 2 P states. 

Table fVIIII is a comparison of the scalar and tensor dipole polarizabilities and hyperpolarizabilities a\, aJ T , 70, 
and 72 for the 2 2 P state of °°Li. In general, our calculations provide significantly more accurate results, particularly, 
for ol^\ A more extensive tabulation of previous results is given in Ref. [27| . 

It should be mentioned that the intermediate configuration of symmetry (pp')P, which contributes to ot\ at the 
level of 0.05%, was not included in the CI-Hylleraas calculation of Pipin and Bishop [42T ]. The relativistic results 
for ai by Wansbeek et al. [l8| and by Johnson et al. [17| are the (2 J + l)-weighted sums between the J = 1/2 and 
J = 3/2 sub-levels. The results of Johnson et al. are close to ours. In contrast, the results of Wansbeek et al. [HI 

deviate significantly from ours; for example, the value of reported by Wansbeek et al. [l8| is a factor of 3.6 larger 
than our calculation. 

The uncertainties in the experimental values for a\ and obtained by Hunter et al. p3| and by Windholz et 
al. [44| are too large to reveal finite nuclear mass and relativistic effects. However, the nonrelativistic calculations of 
Refs7[42j], and the present work, as well as the relativistic calculations of Ref. [l7|, agree with the experiment of 
Windholz et al, though in contrast, as was also observed by Johnson et ai, the values obtained using RKR-QDT by 
Cohen and Themelis [271 ] disagree with the experiment. 

For the ai and 013 of the 2 2 P state, the model potential results of Zhang et al. 16j are, respectively, 4104.9 and 
3.2135 x 10 5 , which are slightly larger than our results, given in Table rVHl 

There are no measurements of 70 and 72 though our calculated values are in almost perfect agreement with, though 
substantially more accurate than, the semi-empirical results of Cohen and Themelis (2?| and the early calculation 
of Themelis and Nicolaides [45[ , obtained by fitting electric-field induced energy shifts calculated using Hartree-Fock 
wave functions. 
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C. Long-range interactions between a 2 2 S atom and a 2 2 P atom 



Table ITXl lists the final values of C3, Cq, and C% for the system 2 2 S-2 2 P between two like-atoms °°Li, 7 Li, and 6 Li 
with all possible symmetries. 

Table 1X1 contains comparisons for C3, Cq, and Cg for the system °°Li(2 2 S)-°°Li(2 2 P) with some published re- 
sults, including the Hylleraas-type calculations of Yan et al. [12j |. the model potential approach of Marinescu and 
Dalgarno [46| . and the semi-empirical model potential of Zhang et al. 16]. Our ab initio results confirm the more 
accurate semi-empirical results of Zhang et al. [1 61 ] , though their results are systematically slightly larger in magnitude 
than ours. It is also evident that the present results have substantially improved the precision of C§. 



Tables IXII and IXII1 respectively, list convergence studies for the dipole polarizabilities ct\ and a \ ' and for the 
hyperpolarizabilities 70, 72, and 74 of 00 Li in 3 2 D. 

Table IXTfll is the summary of all the values for the polarizabilities and the hyperpolarizabilities of °°Li, 7 Li, and 
6 Li in their 3 2 D states. 

Table IXTVl presents comparisons of our calculated scalar and tensor dipole polarizabilities of °°Li with the existing 
theoretical and experimental values for the 3 2 D state. The relative signs and magnitudes (though of limited accuracy) 
were correctly predicted by Themelis and Nicolaides [47j using an empirically modeled dipole moment operator and 
fitting to field-dependent energy shifts. For a±, the percentage difference between the semi-empirical model potential 
result of Magnier and Aubert-Frecon [48[ and ours is about 7%. The values of Ashby and van Wi.jngaarden 49], 
obtained using the semi-empirical Coulomb approximation method and the values of Zhang et al. [16J obtained in 
the framework of a frozen core Hamiltonian with a semi-empirical model potential, are very close to each other, but 
they differ from our results at the levels of, respectively, 0.9% and 0.8%. The experimental results of Ashby et al. [50j 
and the relativistic calculations of Wansbeek et al. [18| are expressed in the LSJ coupling scheme, which may be 
converted into the LS coupling by averaging over the fine structure using a (2 J + l)-weighted sum [lfj]. One can see 
that our result for ol\ disagrees with the experiment at the level of 1%. Similarly, for a\ , the discrepancy is at the 
level of 0.5%. It should be pointed out that, at this level of accuracy, the finite nuclear mass and relativistic effects do 
not account for the discrepancy. The present calculation and that of Wansbeek et al. p~8j are ab initio. Their results, 
which may suffer from an incomplete treatment of correlation effects, differ drastically from our own. 



As by-products in calculating the dipole polarizabilities, we have obtained the oscillator strengths of °°Li, 7 Li, 
and 6 Li for the transitions of 2 2 S - 2 2 P, 2 2 P - 3 2 L>, and 3 2 D - 4 2 F, listed in Table [XV] For the case of °°Li, a 
comparison with some previous results is presented in Table IXVI1 including the multiconfiguration Hartree-Fock of 
Godefroid et al. (l5| , the semi-empirical model potential of Zhang et al. (l6j , and the relativistic many-body approach 
of Johnson [l7| . The nonlinear variational parameters in our calculation are optimized only for the lowest energy 
eigenstate of given symmetry and thus individual oscillator strengths for excited state transitions may not necessarily 
be of better accuracy than other dedicated calculations. Nevertheless, the present calculations have slightly improved 
the previous values [2(| for the 2 2 S — 2 2 P and 2 2 P — 3 2 D transitions, which were obtained variationally using smaller 
Hyllcraas basis sets up to about 3500 terms. The most accurate results for the 2 2 S — 2 2 P oscillator strength was 
obtained by Yan et al. [5l| using variational trial functions that contained the core and valence electron wave functions 
in the zeroth order. [Note that for 7 Li and 6 Li, the present definition for the oscillator strength Eq. f4"3")) differs from 
the one adopted in Refs. 0, [STj] by a factor of 1 + (3 /mo), where too is the nuclear mass.] 

Finally, we have calculated the oscillator strength sum rule [52| S(— 3) = 2/3i, where /?i is given in terms of the 
quantities defined in Sec. [Til m a manner similar to the definition of cti in Eqs. (|3"T)) and (f55|) . as 



D. 3 2 D state: Polarizabilities and hyperpolarizabilities 



E. Some oscillator strengths and sum rules 




(57) 



where 




(58) 
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The results for ft are given in Table IXVIII for °°Li, 7 Li, and 6 Li in the 2 2 S, 2 2 P, and 3 2 D states. The result for 
°°Li in the 2 2 S state is in agreement with the value ft = 1197 given by Pipin and Bishop j37j . 

Another sum rule S(—l) can be related to the interaction potential between an atom and a perfectly conducting 
wall, see Ref. (53[. A highly accurate value obtained [53| using matrix elements from Ref. [54| is 12.144 004 08(24) and 
the present work does not give an improvement. The value from Bishop and Pipin [37] ] is 12.13. 



IV. CONCLUSIONS 

In this paper, the nonrelativistic polarizabilities and hyperpolarizabilities have been calculated ab initio in a unified 
manner for the 2 2 S, 2 2 P, and 3 2 D states of lithium using fully-correlated Hylleraas basis sets. The dispersion 
coefficients for Li(2 2 5)-Li(2 2 S) and Li(2 2 S*)-Li(2 2 P) have also been evaluated. Furthermore, the finite nuclear 
mass effects on these properties have been studied for Li laying the foundation for future work such as investigating 
relativistic and QED effects on the polarizabilities, especially on the dipole polarizability of the ground state of lithium, 
using the approach of Pachucki and Sapirstein [3y| , following the progress on He [36|, [5f| [5(| [57} . 

Our results can also be used as a benchmark for other methods that may be developed in future research. For the 
dipole polarizabilities of lithium in 3 2 D state, an improved measurement would be important in resolving the existing 
discrepancy between the experimental values of Ashby el al. [50| and the present results. 
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TABLE I: Convergence of Qi for the ground state 2 2 S of °°Li, in atomic 
units. 



(N (h Np) ai 

(120, 55) 163.92934 
(256, 138) 164.03473 
(502, 306) 164.06958 
(918, 622) 164.09123 
(1589,1174) 164.10058 
(2625,2091) 164.10695 
(4172,3543) 164.11021 
(6412,5761) 164.11154 
Extrapolated 164.112(1) 



TABLE II: Convergence for the hyperpolarizability of °°Li in the ground 
state 2 2 S, in atomic units. 

~JW,N s ,Np,N D ) T(1,0,1) T(l,2,l) 7(1 

(120,120,55,55) -31221.470 78063.181 177.926 

(256,256,138,138) -31195.797 78088.509 1853.126 

(502,502,306,306) -31177.012 78100.613 2958.636 

(918,918,622,622) -31176.318 78103.456 3044.300 

(1589,1589,1174,1174) -31177.339 78104.691 3019.625 

(2625,2625,2091,2091) -31176.724 78102.359 3004.784 

(4172,4172,3543,3543) -31175.663 78100.447 3018.655 

(6412,6412,5761,5761) -31174.753 78099.257 3038.921 

Extrapolated 3060(40) 
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TABLE III: Values of the polarizabilities ati, a.2, and «3, the hyperpo- 
larizability 70, and the dispersion coefficients Ce, Cs, and C10 for 00 Li, 
7 Li, and 6 Li in their ground states 2 2 S, in atomic units. 



System ai 02 aa 70 Ce Cs C10 

°°Li 164.112(1) 1423.263(3) 39649.29(2) 3060(40) 1393.42(5) 83429(1) 73725(2) x 10 2 

7 Li 164.161(1) 1423.415(5) 39653.72(3) 2820(40) 1394.05(5) 83456(5) 73742(2) xlO 2 

6 Li 164.169(1) 1423.439(4) 39654.46(3) 2780(40) 1394.16(5) 83460(5) 73745(2) x 10 2 



TABLE IV: Comparison of the polarizabilities ai, 0:2, and 03, and the 
dispersion coefficients Ce, Cs, and C10 for t he g round state 2 2 S of 00 Li, 
in atomic units. The results from Refs.[28l]. [29| [3p|. [l8| ]. and fl7] ] were 
calculated using relativistic methods. 



Reference 



f-U 



02 



Q'3 



c e 



Cs 



C10 



Yan et al. [12] (1996) 
Komasa [26] (2001) 
Derevianko et al. [28] (2001) 
Porsev and Derevianko [29] (2003) 
Chen and Wang [25] (2004) 
Cohen and Themelis [27] (2005) 
Zhang et al. [16] (2007) 
Sahoo [30] (2007) 
Wansbeek al. [18] (2008) 
Johnson et al. [17] (2008) 
This work 



39650.49(8) 1393.39(16) 8.34258(42) x 10 4 7.3721(1) x 10 b 
1389(2) 



164.111(2) 
164.11171 
164.0(1) 



164.14 

164.21 

162.48(56) 

162.87 

164.084 

164.112(1) 



Molof et al. [32] (experiment) (1974) 164.0(3.4) 
Miffre et al. [31] (experiment) (2006) 164.2(1.1) 



1423.266(5) 
1423.282 

1424(4) 
1423.48(17) 
1423.3 
1424.4 

1421.37(3.51) 
1420 
1422.73 
1423.263(3) 



3.957xl0 4 
39650.96(94) 

39680 



39624.2 
39649.29(2) 



1394.6 

1396(6) 

1393.42(5) 



8.34(4) xlO 4 7.35 xlO 6 



8.3515xl0 4 



8.360xl0 4 



7.3811 xl0 b 



8.3429(1) xlO 4 7.3725(2) xlO 6 



TABLE V: Values of the hyperpolarizability 70 for °°Li in the ground 
state 2 2 S, in atomic units. 



Reference 70 

Pipin and Bishop [37] (1992) 3000 

Kassimi and Thakkar [38, 39] (1994) 2900(90) 

Laughlin [40] (1995) 3930 

Jaszuhski and Rizzo [41] (1996) 3450 

Cohen and Themelis [27] (2005) 3390 

Present 3060(40) 



TABLE VI: Convergence of ai and for °°Li in 2 2 P state, in atomic 
units. 



(N ,Ns,N (pp , )P ,N D ) ai (S) ai(( PP ')P) ai(D) 



a 1 



TJTT 



(55,120,36,55) -17.181258 0.063053 142.611225 125.493021 2.951662 

(138,256,108,138) -16.258204 0.067825 142.675018 126.484639 2.024615 

(306,502,264,306) -15.950171 0.069135 142.718553 126.837517 1.712883 

(622,918,568,622) -15.871922 0.069501 142.729209 126.926788 1.633752 

(1174,1589,1106,1174) -15.866615 0.069591 142.734228 126.937203 1.627988 

(2091,2625,2002,2091) -15.861376 0.069641 142.735430 126.943695 1.622654 

(3543,4172,3413,3543) -15.860744 0.069646 142.736141 126.945043 1.621953 

(5761,6412,3413,5761) -15.860549 0.069648 142.736441 126.945540 1.621729 

Extrapolated 126.9458(3) 1.6214(3) 
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TABLE VII: Values of the polarizabilities ai, , a 2 , and a 3 , and the 
hyperpolarizabilities 70 and 72 for °°Li, 7 Li, and 6 Li in 2 2 P state, in 
atomic units. 



System 


<*i a\ a 2 03 7o 


72 


°°Li 
7 Li 
6 Li 


126.9458(3) 1.6214(3) 4103.165(5) 321138(4) 1.00170(9) xlO 7 
126.9472(5) 1.6351(2) 4102.893(4) 321102(5) 1.00201(9) x 10 7 
126.9474(5) 1.6373(3) 4102.848(4) 321096(5) 1.00206(9) x 10 7 


-6.2234(8) xl0 B 
-6.2252(8) xlO 6 
-6.2255(8) xlO 6 



TABLE VIII: Comparison of the scalar and tensor dipole polarizabilities 
and hyperpolarizabilities for °°Li in 2 2 P state, in atomic units. The 
results from [l8l ] and [l7| ] are relativistic. 



Reference 


ai 


af> 


7o 


72 


Themelis and Nicolaides [45] (1992) 


135.7 


0.13 


1.10x10' 


-6.970 xl0 b 


Pipin and Bishop [42] (1993) 


126.844 


1.605 






Cohen and Themelis [27] (2005) 


126.4 


1.73 


1.002xl0 7 


-6.21 xlO 6 


Zhang et al. [16] (2007) 


126.95 


1.6627 






Wansbeek et al. [18] (2008) 


125.20 


5.95 






Johnson et al. [17] (2008) 


126.990 


1.59 






This work 


126.9458(3) 1.6214(3) 1.00170(9) xlO 7 


-6.2234(8) xlO 6 


Hunter et al. [43] (1991) (experiment) 


126.8(3.4) 









Windholz et al. [44] (1992) (experiment) 126.87(36) 1.64(4) 



TABLE IX: Values of C 3 , C 6 , and C 8 for 2 2 S-2 2 P of ' 
in atomic units. 



3 Li, 7 Li, and "Li, 



System C 3 (M 2 = 0) C 3 (M 2 = ±1) C 6 (M 2 = 0) C a (M 2 = ±1) C 8 {M 2 = 0) C S {M 2 = ±1) 



°°Li 
7 Li 
6 Li 



OOT ■ 

Li 

7 Li 
6 Li 



-11.000221(2) 5.500111(1) 

-11.001853(2) 5.500926(1) 

-11.002125(2) 5.501062(1) 

11.000221(2) -5.500111(1) 

11.001853(2) -5.500926(1) 

11.002125(2) -5.501062(1) 



H = -l 
2075.40(3) 
2076.08(7) 
2076.19(7) 

= +l 
2075.40(3) 
2076.08(7) 
2076.19(7) 



1406.68(3) 990895(5) 48564.8(5) 

1407.15(5) 991075(6) 48566.4(2) 

1407.20(2) 991104(5) 48566.9(4) 

1406.68(3) 274079(2) 103044(2) 

1407.15(5) 274128(5) 103052(1) 

1407.20(2) 274137(6) 103053(1) 



TABLE X: Comparison of C 3 , C 6 , and C 8 for the system °°Li(2 2 5')- 
°°Li(2 2 P), in atomic units. 



Reference f3 C 3 (M 2 = 0) C 3 (M 2 = ±1) C 6 (M 2 = 0) C 6 (M 2 = ±1) C 8 (M 2 = 0) C 8 (Af 2 = ±1) 



Marinescu et al. [46] (1995) 


-1 
+1 


-11.01 5.503 2066 


1401 


9.880 x 10 b 
2.705 x 10 5 


4.756 x 10 4 
1.021 x 10 5 


Yan et al. [12] (1996) 


-1 


-11.000226(15) 5.5001133(74) 2075.05(5) 


1406.08(5) 


9.9202 x 10 5 
2.7431 x 10 5 




Zhang et al. [16J (2007) 


-1 
+1 


-11.008 5.5041 2076.3 


1407.4 


4.8629 x 10 4 
1.0316 x 10 5 


This work 


-1 
+ 1 


-11.000221(2) 5.500111(1) 2075.40(3) 


1406.68(3) 


9.90895(5) x 10 5 
2.74079(2) x 10 5 


4.85648(5) x 10 4 
1.03044(2) x 10 5 



TABLE XI: Convergence for otx and a^ T ' of °°Li in 3 2 D state, in atomic 
units. 



(N ,N P ,N (pd)D ,N F ) 



MP) 



Ql(( P d)D) Qi(F) 



Ql 



AT) 



(138,138,126,132) 
(306,306,322,302) 
(622,622,714,636) 
(1174,1174,1428,1248) 
(2091,2091,2640,2307) 
(3543,3543,4587,4051) 
(5761,5761,4587,6806) 



48857.791600 
-16933.479208 
-16850.080546 
-16845.626886 
-16845.378158 
-16845.342790 
-16845.342870 



0.059150 
0.062827 
0.063886 
0.064157 
0.064226 
0.064242 
0.064243 



1916.643488 
1916.940725 
1917.021571 
1917.026639 
1917.035892 
1917.040003 
1917.043303 



-16941.088962 
-15016.475656 
-14932.995089 
-14928.536090 
-14928.278040 
-14928.238545 
-14928.235324 



18310.238324 
16385.844685 
16302.423983 
16297.969146 
16297.717843 
16297.681317 
16297.680455 



TABLE XII: Convergence for 70, 72, and 74 of °°Li in 3 2 D state, in 
atomic units. 



(N ,Ns,Np,N (pp , )P ,N D ,N (pd)D ,N F ,N (pf)F ,N G ) 10^70 10^72 10~ lu 74 



(138, 256, 138, 108, 138, 126, 132, 126, 139) 
(306, 502, 306, 264, 306, 322, 302, 322, 330) 
(622, 918, 622, 568, 622, 714, 636, 714, 720) 
(1174,1589,1174,1106,1174,1428,1248,1428,1458) 
(2091,2625,2091,2002,2091,2640,2307,2640,2769) 
(3543,4172,3543,3413,3543,4587,4051,4587,4975) 



2.335078929997 -2.435345979367 10.1097221264 

1.673627590440 -1.744707390466 7.1916589358 

1.648569007586 -1.718605330888 7.0875021256 

1.647221992572 -1.717197712093 7.0814607892 

1.647149117913 -1.717122212016 7.0812009263 

1.647145129442 -1.717120442563 7.0814240701 



TABLE XIII: Values of the polarizabilities at, a\ ,0.2, and 0:3, and the 
hyperpolarizabilities 70, 72, and 74 for °°Li, 7 Li, and 6 Li in 3 2 D state, 
in atomic units. 



System 


(T) 

Qi a\ a.2 


0:3 


7o 


72 


74 


00 Li 
7 Li 
6 Li 


-14928.230(5) 16297.675(5) 158060(10) 
-14921.330(4) 16291.094(5) 158070(10) 
-14920.180(6) 16290.000(5) 158070(10) 


-1.340902(3) x 10 s 
-1.339746(5) x 10 8 
-1.339554(5) x 10 s 


1.647140(5) x 10 12 
1.644875(5) x 10 12 
1.644500(4) x 10 12 


-1.717115(5) x 10 12 
-1.714740(5) x 10 12 
-1.714345(5) x 10 12 


7.0814(2) x 10 1U 
7.0700(3) x 10 10 
7.0680(5) x 10 10 



TABLE XIV: Comparison of the scalar and tensor polarizabilities and 
hyperpolarizabilities for °°Li in 3 2 D state, in atomic units. The results 
from T8J are relativistic. 



Reference 


10" 4 «i 


io-M T) 


10~ 5 «2 


10" 8 a 3 10" 12 7o 


10" 12 72 10" 10 74 


Themelis and Nicolaides [47] (1995) 


-2.0468 


2.1944 




4.56 


-3.97 19.5 


Magnier and Aubert-Frecon [48] (2002) 


-1.3950 


1.5324 








Ashby and Wijngaarden [49] (2003) 


-1.507 


1.642 








Zhang et al. [16] (2007) 


-1.5044 


1.6414 


1.5786 


-1.3548 




Wansbeek et al. [18] (2008) 


-1.986 


2.090 








This work 


-1.4928230(5) 1.6297675(5) 1.58060(10) 


-1.340902(3) 1.647140(5) 


-1.717115(5) 7.0814(2) 


Ashby et al. [50] (2003) (experiment) 


-1.513(4) 


1.643(6) 









TABLE XV: Values of oscillator strengths of °°Li, 7 Li, and 6 Li. 



System 2 2 g-2 2 P 2 2 P-3 2 D 3 2 £>-4 2 F~ 

°°Li 0.7469563(5) 0.6385685(5) 1.0153771(5) 

7 Li 0.7469614(4) 0.6385835(5) 1.0154562(5) 

6 Li 0.7469623(4) 0.6385858(4) 1.0154695(5) 



TABLE XVI: Comparison of oscillator strengths of °°Li. 
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Reference 2 2 S - 2 2 P 2 2 P - 3 2 D 3 ' 2 D - 4 2 F 

Yan and Drake [20] (1995) 0.7469572(10) 0.6385705(30) 

Yan et al. [51] (1998) 0.7469569396(98) 

Godefroid et al. [15] (2001) 0.74690 0.63853 

Zhang et al. [16] (2007) 0.7475 0.6388 1.0153 

Johnson et al. [17] (2008) 0.746944 0.638615 1.015637 

This work 0.7469563(5) 0.6385685(5) 1.0153771(5) 



TABLE XVII: Values of ft of °°Li, 7 Li, and 6 Li in their 2 2 S, 2 2 P and 
3 2 D states, in atomic units. 



System ft (2^) ft(2"P) ft(3~^D)~ 

°°Li 1196.9696(2) 1614.68(2) 5202428(1) 

7 Li 1197.4886(2) 1614.99(2) 5197637(1) 

6 Li 1197.5750(2) 1615.04(2) 5196841(1) 



APPENDIX A: STARK EFFECT 



The Hamiltonian for an atom in a uniform electric field £ = £ z is given by 

H = H Q + H' = H -e-P, (Al) 
where Ho is the unperturbed Hamiltonian and P is the electric dipole moment of the atom: 

p = • ( A2 ) 

i 

In the above, qi is the charge of ith particle, R; is its position vector relative to a laboratory frame, and the summation 
is over all charged particles inside the atom, including the nucleus. Under the perturbation H', the energy eigenvalue 
and cigcnfunction of H can be written in the form 

E = E + AE X + AE 2 + AE 3 + AE t + • • • , (A3) 
# = + *i + *2 • • • , (A4) 



where 



H *o = S *o (A5) 



is the zero-order equation, and AEi and ^ are the corresponding ith-order corrections. According to the perturbation 
theory, the energy corrections can be expressed as 

AE 1 = (* \H'\* ) , (A6) 

AE 2 = (* |#'|*i), (A7) 

AE 3 = -AEi (A8) 

AE 4 = (#i|ff'|#2)-A£ 2 (#i|#i)-A.Ei(#i|#2). (A9) 

If the state of interest * has a fixed parity, as in the case of this work, AEi — due to the parity selection rule. *i 
and $2 can be expanded in terms of their spectral representations: 



m\ X- (n\H'\k)(k\H'\0) 

'* 2> = ^(E -E n )(E -E k )\ n) > (AU) 

nk 
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where |0) = |\&o) and {E n , \n)} is a complete set of Hq, including the continuum. Inserting Eq. (|A10[) into Eq. (|A8 
yields A£" 3 — 0, also due to the parity consideration. Thus, 



AE 2 = 

n 

£ 



|H>)(4tf]0) 
Eq — E n 



AE. 



kmn 



(Eq — E m )(E ~ E n )(E Q — Ek) 



S(n,0) 



1 



(E - E m )(E - E k f 



x(0\H'\m}(m\H'\n)(n\H'\k)(k\H'\0) 



(A12) 



(A13) 



1. AE 2 

Let us first consider the operator H'\n){n\H' in Eq. (|A12|) . Using spherical tensor operator technique, we have the 
following decomposition 

H'\n){n\H' = ^(-1) A '+«[P« ® A„P«] <f>[E (1) ® S^] (K q ] , (A14) 

Kq 

where A = | tt,} (77, | . Since £ is along the z-axis, only the q = component survives in the above equation. Thus, 

H'\n)(n\H' = £(-1)^2^+1 £ Q * f) P«A n P«[^) ® > . (A15) 

if 9192 

After substituting the Eq (|A15[) into Eq (|A12|) . one has 

a Ej = £(-i)*v5FTT£ 0, i, 0) . (am) 

nK 9192 ^ ' On 

To be specific, let us write out explicitly the angular momentum quantum numbers in the initial and intermediate 
states: 

|0) = \n LM), (A17) 
\n) = \nL a M a ), (A18) 

where no and n are the corresponding principal quantum numbers. Then, the summation over n in Eq (|A16|) actually 
means the summation over {n, L a ,M Q }. By applying the Wigner-Eckart theorem for the irreducible tensor operator 
P^p , one can recast Eq (|A16p into the following form 



AE 2 = VM^vWTT ^ 11 ^ 1 ^ 

En (E) — E n (L a ) 



-iL„K 



(A19) 



where E no (L) and E n (L a ) stand for E and E n respectively in Eq. (|A16[) . and 



A _ E _^ ; _^ ; . ; x « . (A20) 

The quantity A can be simplified by using the standard graphical method of dealing with angular momentum [58| , 
provided K is an integer, and note that 

[*« eH K > = VSr+T £ ( £ * f) = V2ZT1 (J J f) £ 2 ■ (A22) 



Therefore, the second-order correction can be written 

(noL\\pW\\nL a )(nL a \\pW\\n L) 



AE 2 = -E 2 J2 ' 



7lL a 



E n (L a ) — E no (L) 



B-^)(-')'-"(:^)(- L M^^){iif„ 



Since 



x 

K 



\L-M I L L 0\ / i \2L ' 



( 1} \-M M oj - ( X) V2LTT' 

/^t-tf/ L L 2\ 2i 3Af 2 -LQL + l) 

1 j M V ( ' v /(2i + 3)(L + l)(2L + l)L(2L-l)' " 



and also [5S 



p i 1} = yYE^^CRi), 

(nL ||pW||noi) = (-1) L - L ° (n L||P«|K a >* 
the second-order energy correction can finally be expressed in the form 

AE 2 = ~{a 1 +af ) g 2 {L,M)). 

In the above, g 2 (L, M) is the only M-dependent part: 

0, if L = 0, \ 

g 2 (L,M) = { 3M 2 -L(L+1) 



L(2L - 1) 



- , otherwise 



(T) 

and a\ and a\ are, respectively, the scalar and tensor dipole polarizabilities 

a[ T ) = ^ L a )a\ (L a ) , 



where 



with Ti = £\ £?i-Ri^io(R»)) and 



11 qJ 9(2Z + 1) 4- -£ no (L) 



W ,L a ) = ( _ 1) ^,/30(2L +1 )X(2X-l)ri 1 2 



y (2i + 3)(i + l) \i £ L a 
2. AE 4 

According to Eq. (|A15[) , the fourth-order energy correction of Eq. (|A13|) can be written as 

AE, = £t(fc,m,n) £ {-ir^iK^K^ £ 
kmn KtK 2 qiqiq-s 





1 




1 


1 






'72 




\93 


?4 
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where the notation (a, b) = (2a + 1)(2& + 1), 

^ ^ (Po — E m )(E — E n )(E — E k ) ^ ' \e — E m )(E — E k ) 2 ' ^ ^ 

and A m = |m)(77i|, etc.. Writing out the all angular momentum quantum numbers explicitly 

|0) = \n LM), (A36) 

\m) = \mL a M a ), (A37) 

|n) = \nL b M b ), (A38) 

|fc) = \kL c Mc), (A39) 

and applying the Wigner-Eckart theorem, we have 

AP 4 = E E t(k,m,n) ^ (-l)^^,^) 1 / 2 ^ 1 ' ^W]^'^ 1 ) ^W]f s) 

kmn L a LbL c K\Ki 

x (rioillpWllmL^^iallPWllri^^XfellpWll^^LcllpWllnoL)^, (A40) 
where P (which contains all the angular coefficients) is 

v L - M+L a - M a +L b - M b +L C - M c 



b= E E (-!) J 

9192 93 94 M a M b M a 

L a 1 LiA A L h I L c \( L c I L 
-M a qi M b J y~M b q 3 M C J \-M c q 4 M 

Using the graphical method, B can be simplified into 





i 


K 2 \ 


{ L 


1 




Hi 


?4 


o ) 


\-M 


<h 


Ma) 



(A41) 



L L b L a \L L b L c ^ MO 0/ I -M M 0/ I £ £ 
provided K\ and -ft^ are integers. From Eqs. (|A22[) , (I A26|) . and (|A27|) . one can further write Eq. (|A40|) in the form 

AP 4 = E nL a ,L b ,L c )J2(-l)^(_ L M L M J) ^(A,^,^ 

(A43) 



-M M 0, 

v ~ 7 L a L b L c A 1 7 K"iK 2 

X 2 A\ f 1 1 Kx\ f 1 1 if 2 \ / ^2 A"i A 
\ L L b L a [ 1 L L b L c f 1 L L L b 





1 






i 




G 







G 





o J 



where 



Tl7 T j"| = V ("-o ^ 1 1 1 1 ra£ Q ) (raL a 1 1 ^1 1 1 (nlj 1 1 J\ 1 1 kL c ) (kL c | | Ty 1 1 n L) 

" b ' " fc™ [P fc (L c )-P„o(i)][^ m (ia)-P, ln (i)][P„(i b )-P„oW] 

" 5 ^ L){ - 1) 1^ [Em (L a ) - E no (L)]\ [E k (L c) - E no (L W ■ (M4) 



According to the property of 3j symbol, the possible values for A are 0, 2, and 4. Also, 
(-!)> 



,l-m L L 4 
y-M M 

2(-l) 2L [3(5M 2 — L 2 — 2L)(5M 2 + 1 - £ 2 ) - 10Af 2 (4M 2 - 1)] 
a/(2Z + 5) (2L + 4)(2L + 3)(2L + 2){2L + 1)(2L)(2L - 1)(2L - 2){2L - 3) ' ~ ' 

Together with Eqs. (|A24p and (|A25|) . the fourth-order correction can finally be expressed in the form 

£ A 

A£ 4 = (70+72 g 2 (L,M))+ 74g4 (L,M)), (A46) 
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where g2(L,M) is defined in Eq. (|A29[) . and §i{L,M) is given by 

gJL, M) = { 3(5M 2 — L 2 — 2L)(5M 2 + 1 - L 2 ) - 10M 2 (4M 2 - 1) , . (A47) 
' otherwise 



0, if £ < |, 



L(2L- l)(2L-2)(2L-3) 



In Eq. (|A46[) . 70 is the scalar hyperpolarizability, and 72 and 74 are the tensor hyperpolarizabilities, which can be 
written as 



;n : (-l) 2L l2 ^ 2 Af ^= V g a (L,L a ,L b ,L c )T(L a ,L b ,L c ), ( A 4s ) 

" ^ L„L h L a 



7 2 = (-1)"^ J (arw f+'i) 1 ^ + 1) S (A49) 

V L a LijL c 



where 



= ( _ 1)2i 128^ / y- l)(£-l)(2L- 3) ^ Q^(L, L a , Lj,, L c )T{L a , Lb, L c ), (A50) 

' 1 ; 3 V (2L + 5 L + 2 2L + 3 (L + 1 (2L + 1 r , r 

Ga(L, L a , L b , L c ) = J] (A.Jfx.Jfa) (j J ^) (j J A q 2 ) ^ q) {i L 6 £} 



1 1 isr 2 l fi^ 2 ifi a 



iv Z/5 Z/ c f 1 Z/ iv Li, 



(A51) 



